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A generalized Gross-Pitaevskii equation adapted to the U( 5) D SO( 5) D SO( 3) symmetry has 
been derived and solved for the spin-2 condensates. The spin-textile and the degeneracy of the 
ground state (g.s.) together with the factors affecting the stability of the g.s., such as the gap and 
the level density in the neighborhood of the g.s., have been studied. Based on a rigorous treatment of 
the spin-degrees of freedom, the spin-textiles can be understood in a IV-body language. In addition 
to the ferro-, polar, and cyclic phases, the g.s. might in a mixture of them when 0 < M < 2N (M 
is the total magnetization). The great difference in the stability and degeneracy of the g.s. caused 
by varying p (which marks the features of the interaction) and M is notable. Since the root mean 
square radius J? rms is an observable, efforts have been made to derive a set of formulae to relate 
Rims and N, ^(frequency of the trap), and p. These formulae provide a way to check the theories 
with experimental data. 
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I. INTRODUCTION 

Since the success of trapping cold atoms via optical 
trapjlj the study of condensates of atoms with nonzero 
spin is a hot topic. There are numerous literatures 
dedicated to the study of the ground state (g.s.). For 
spin-2 condensates, the g.s. was found to have three 
phases, namely, the ferro-phase, polar phase and cyclic 
phase, depending on the features of the interaction. pj-Q 
On the other hand, it has been revealed that the total 
Hamiltonian of the spin-2 systems is associated with the 
U( 5) D 50(5) D SO( 3) symmetry. (H, @ Since there are 
elegant theory dealing with this symmetry, it is worthy 
to derive a generalized Gross-Pitaevskii (GP) equation 
adapted to the 1/(5) D 50(5) D 50(3) symmetry. Since 
the solutions of this equation could provide the details of 
the low-lying states, we believe that the understanding 
to the spin-2 condensates could be thereby enriched. 

We consider the case that the magnetic field B = 0. 
The symmetry-adapted GP equation is derived firstly. 
Then, by solving the GP equation, both the spin-textiles 
and the spatial wave functions of the low-lying states 
can be known. The emphasis is placed on the g.s.. Since 
the stability of the g.s. is affected by its environment, 
in addition to the g.s. itself, special effort is made to 
study the environment, namely, the width of the gap (the 
energy difference between the first excited state and the 
g.s.) and the level-density in the neighborhood of the 
g.s.. This is a topic less studied before. Besides, the 
degeneracy of the g.s. is an important feature which is 
also studied in this paper. The total magnetization M 
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is considered as being conserved, which depends on how 
the system is experimentally prepared. The influence of 
M on the spin-textiles and the degeneracy of the g.s. is 
found to be great. 

After we have obtained the spatial wave functions, 
physical quantities related to the spatial degrees of free¬ 
dom can be known. In particular, we have calculated 
the root mean square radius R rms . Since R rms is an ob¬ 
servable, theoretical results and experimental data can be 
compared with each other. This is a notable way to check 
the theory. The GP equation will be firstly solved analyt¬ 
ically under the Thomas-Fermi approximation (TFA). A 
number of formulae relating the physical quantities have 
been thereby obtained, and the underlying physics can 
be understood in an analytical way. Then, strict numeri¬ 
cal calculations are performed to check the validity of the 
TFA. 


II. INHERENT SYMMETRY AND THE 
GROSS-PITAEVSKII EQUATION 

A condensate with N spin-2 neutral atoms trapped by 
an isotropic harmonic potential f mu 2 r 2 is considered, 
in which the spin-orbit coupling is assumed to be weak. 
It is further assumed that the trap is not so weak (say, 
u> > 100s _1 ) and the spin-dependent interaction is weak 
so that the spin-modes are much lower than and sepa¬ 
rated from the spatially excited modes. This is the basic 
assumption of this paper. In other words, the following 
results hold only for the systems adapted to this assump¬ 
tion. In this case a group of excited states distinct in 
their spin-modes would emerge in the neighborhood of 
the g.s.. In each of these states (including the g.s.) all 
the particles have a common spatial wave function which 



2 


is most advantageous to binding and adapted to a spe¬ 
cific spin-textile. Thus one of these states with a specific 
spin-mode 7 can be in general written as 

N 

^7 = 07( r ./)©7, (1) 

i=l 

where 0 7 is the common spatial wave function and © 7 is 
a total spin-state which is now unknown. If the particles 
did not fall into the same state but have some of them 
getting excited, the total energy would be thereby higher. 
The group {’P 7 } totally form a band which is called the 
ground band. The states not in the ground band will 
contain various spatially excited modes and will not be 
studied in this paper. _ 

When hco and A = yfh/(mu>) are used as units for 
energy and length (where m is the mass and is given at 
the one of 87 Rb), the total Hamiltonian is 

# = £*« + £ Vy, (2) 

i i<j 

where h{i) = -|V? + \r\. Vij = 6fa - rj) £ s 3 ^'- 
P s j = X(*) IS the 

spin-state of the i-tli particle, and the two spins of i and 
j are coupled to the combined spin s and its z-component 
m s , s = 0, 2, or 4. Obviously, P l J is the projector onto 
the s-channel. g s is the strength related to the scattering 
length of the s-channel. 

We introduce an operator defined in the total spin- 
space as V e . ]T) s g s Ps 3 ■ It has been proved that 

V is a linear combination of a set of Casimir opera¬ 
tors belonging to a^chain of nested algebra as U( 5) D 
50(5) d 50(3). H, @] Consequently, the eigenstates and 
the eigenenergies of V can be exactly known as 

V 0 7 = E vS e-y. (3) 

Now, 7 represents a group of quantum numbers v, ntri, 
5, and M, where 5 and M are the total spin and its Z- 
component, v and ntri will be explained below, and © 7 = 
©untriSM is the related eigenstate. For convenience, M > 
0 is assumed. 

When two spins are coupled to zero, they form a sin¬ 
glet pair 0 p air(ij) = (x(*Mi))o- When three spins are 
coupled to zero, they form a triplet and can be approx¬ 
imately denoted as 0 + r i ~ (xxx)o (this notation is ex¬ 
act when N —> 00 ). [7] In the triplet every pair of spins 
are coupled to 2. It turns out that the singlet pair and 
triplet are basic building blocks. An eigenstate may con¬ 
tain 7 ip a i r singlet pairs and ntri triplets. The quantum 
number v = N — 2n pa i r is named the seniority, which 
is the number of particles not in the singlet pairs. The 
number of particles neither in the pairs nor the triplets 
is IVfree = N — 2 n pa ; r — 3n tr i, the total spin 5 is con¬ 
tributed by them. It has been proved that, when n pa i r , 
ntri, 5 and M are given, the totally-symmetric eigen¬ 
state ©tmtriSM is unique. The four quantum numbers 


are bound by the following conditions: (i) Nf lee > 0, (ii) 
-Nfree < 5 < 2Nf iee , (iii) 5 = 27Vf r ee — 1 is not allowed, 
and of course (iv) \M\ < 5.|1, 0 

Let us define g 0 4 = 50 - 34 , 324 = 32 - 34 , and g {024 ) = 
5(30 + 32 + 34 )- Then, the eigenenergy associated with 

©7 is @ 

E v g = d\N + a 2 N(N + 4) + g v v(v -I- 3) -I- gsS (5 +1), (4) 

where g v = (-7 304 + 103 2 4 )/ 70 , g s = (- 324 )/14, ai = 
Tlbo4- ^324 - §3(024), and a 2 = - ^ 304 - ^324 + ^3(024) ■ 
Note that E v s does not depend on ntri implying that the 
levels might be degenerate as shown below. 3 

Inserting Eq. O into the Schrodinger equation, we have 

N N 

<n - e i \n ^(r,)©7)=0, (5) 

3 =2 . 7=1 

where the integration covers all the degrees of freedom 
except ri. From this equation and making use of Eq.(3), 
it is straight forward to obtain the symmetry-adapted 
GP equation for the normalized <^ 7 (ri) as 

[h + —E v s(f>j(fi'y — e 7 ]</> 7 = 0, (6) 

where e 7 is the chemical potential. After solving Eq.( 6 ) 
we can obtain the total energy E 1 = Ne 1 — Vtot, 7 , where 
Vtot , 7 = / |07 | 4 rfr E v s is the total interaction energy. 
The whole spectrum of the ground band can be obtained 
from the set {E 7 }. Obviously, this equation is a general¬ 
ization of the one for spin-1 condensates given in ref. [9) 
It is also a more accurate version for the one given as 
Eq.(74) in ref.Q, in which the spin-dependent interac¬ 
tion has been neglected. 

III. THOMAS-FERMI APPROXIMATION 

Since N is usually large, Eq.(6) can be solved by using 
the TFA. Neglecting the kinetic term in Eq.(6), when Hu> 
and A are used as units, the normalized TFA solution is 

07 = V^ 1 “ ( r Ao) 2 and e 7 = rg/2, where 

r ° = < 2 ^> I/S ’ m 

is the TF-radius. From this solution we can obtain the 
root mean square radius i? r ms = ( r 2 ) 1 / 2 = ^/3/7r 0 , and 
the total energy 

£7 = y(^J) 1/5 fV 3/5 (45) 2/5 = InrL, (8) 

Thus, the total energy is directly related to the size (mea¬ 
sured by R T uis)- 

When N is large, we know from Eq.(4) that E v g is 
proportional to N 2 . Accordingly, i? rms oc A 1 / 5 and E 7 oc 
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FIG. 1: (color on line) The ground state energy E gs (in fico) 
versus ip. N = 12000 and ui = 300s -1 are assumed (the same 
in the following figures except specified). The total magneti¬ 
zation is given at five cases as M = 0, 1000, TV, 2N — 1000, 
and 2N, respectively, for the curves ”1” to ”5”. There are 
three regions of i p separated by ip a = tan -1 ( Ji^n+m/ 2 - 2 ) ) 
and ip b = tan -1 (^y^) marked by the two vertical dotted lines. 
In each region the g.s. has its specific (v, S, d) marked in the 
region, where d is the degeneracy. In region I, E gs is M- 
independent. In I and II, the g.s. is not degenerate (d = 1). 
In III the degeneracy depends on M seriously. For the curves 
”1” to ”5”, d = 1, 167, 2001, 167, and 1, respectively. 

TV 7 / 5 . Thus the size increases with TV very slowly, but Ej 
increases with TV faster than linearly. 

On the other hand, it is noted that, in the units of Huj 
and A, all the strengths g s oc y/Q, where Q is the magni¬ 
tude of u>, i.e., u> = ws _1 . Therefore E v s oc yfcj. When 
w-independent units are used, the size is described by 
i? rms A oc (w) -2/5 , i.e., a stronger trap leads to a smaller 
size. Similarly, one can prove that the chemical poten¬ 
tial e 1 hu oc (w) 6 / 5 , and the total energy E^hcu oc (w) 6 / 5 . 
Thus the total energy will increase with oj a little faster 
than linearly. 


IV. GROUND STATE AND ITS PHASE 

It is obvious from Eq.(4) that which pair of v and 5 
is more advantageous to binding depends on the com¬ 
petition between g v and gs, which depend on go 4 and 
< 724 , and they are considered as variable. Therefore, we 
introduce a variable parameter tp to manifest the com¬ 
petition. In the units of Huj and A, the strengths are 
assumed as <704 = ^ 2 cos ip, 324 = u{Cj) 1 / 2 s\iup, 

ff( 024 ) = 2.5u(d)) 1 / 2 , and u = 10 -3 (For a comparison, 
the realistic strengths for 87 Rb are (724 = — 0.28it(d>) 1 / 2 , 
go 4 = —0.4Tu(cD) 1 / 2 , and #( 024 ) = 2.39u(w) 1 / 2 ). In fact, 
it has been shown previously that the phase of the g.s. 
depends on ip. M 

Under the TFA, the ground state energy E gs is just 
the lowest E 1 given by Eqs.( 8 ) and (4). Although E 1 


does not directly depend on M, the condition 5 > \M\ 
restricts the possible choice of 5, and E gs is thereby af¬ 
fected. E gs versus p with M given at five values are 
shown in Fig.l. There are two critical points ip a = 
tan " 1 ( _ 7 ff^/ 2 ,+ 3) 2 ) ) and ip b = tan-^f^j.Q When ip 
goes through either ip a or ip b a phase transition will occur. 
This phenomenon was found as early as in 2000 by using 
the mean field theory (in which ip a = tan - ^^^,) and the 
same ip b were found .[2 The present more accurate result 
will tend to the old one when N is large). 

When ip is in (0, ip a ), i.e., in region-I, the g.s. does not 
depend on M and will have (v, 5, d) = ( N , 2 N, 1), where 
d is the degeneracy of the state. Since all the spins align 
along a common direction in this g.s., it is in the ferro- 
phase. @-[ 1 ] 

When ip is in (ip a ,<p b ), i.e., in region-II, the g.s. de¬ 
pends on M and will have (v, 5, d) = (M/2, M, 1). Since 
v = M/2 and 5 = M, all the M/2 unpaired spins must 
align along a common direction so that the total spin 
can be maximized (i.e., 5 can be equal M). Thus each 
g.s. in region-II is a mixture of a group of aligned spins 
together with the (N — Mj 2)/2 singlet pairs, namely, a 
mixture of ferro-phase and polar phase. When M = 0, 
(■v , S, d) = (0,0,1), thus it is in pure polar phase (we have 
assumed that N is even, otherwise 5 = 2). When M is 
larger, more spins are in the ferro-phase, and accordingly 
Eg a gets higher as shown in the figure. When M = 2 TV, 
(■v , 5, d) = (TV, 2 N, 1), and the g.s. is in pure ferro-phase. 

When ip is in (ip b , 27t), i.e., in region-III, the g.s. de¬ 
pends also on M and has ( v,S,d ) = ( N,M,d ). Note 
that these g.s. have v = N, therefore they do not con¬ 
tain any singlet pairs, but the spins in triplets are al¬ 
lowed. It is recalled that, in regions I and II, the g.s. has 
5 = 2v. It implies that all the unpaired spins must be 
aligned, and therefore there is no room for the triplets. 
Whereas in III, the number of triplets nto can have dif¬ 
ferent choices in a specific domain under the constraints 
of symmetry given in points (i) to (iv) in the paragraph 
above Eq.(4) (while those spins not in the triplets are 
coupled to 5 = M). This leads to the degeneracy of the 
g.s. [8} When M = 0 the g.s. has 5 = 0 and therefore 
iVf r ee = 0 (point (ii)). Thus, all the particles are in the 
triplets (i.e., in a pure cyclic phase), and henceforth d = 1 
(we have assumed that N is a multiple of 3 to simplify 
the discussion). When M = 2 N, the g.s. has 5 = 2 N 
and therefore N{ lee = TV. Thus, neither the pairs nor the 
triplets emerge (i.e., in a pure ferro-phase), and hence¬ 
forth d = 1 also. When M is neither 0 nor 2TV, the g.s. is 
in general degenerate. In particular, when M = TV, TVf ree 
can be ranged from M/2 to M (each step is 3). In this 
case the degeneracy d is maximized (the curve ”3” in III 
has d = 2001). The appearance of highly degenerate g.s. 
in the region III is a notable feature. 

It is reminded that the polar phase is composed of 
s = 0 pairs. Therefore, when go is less positive than <72 
and 5 ( 4 , the g.s. will prefer this phase. The cyclic phase 
is composed of s = 0 triplets in which every pair of spins 
must be coupled to s = 2. Therefore, when <72 is less 
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positive, the g.s. will prefer this phase. While the ferro- 
phase is composed of aligned spins, in which every two 
must be coupled to s = 4. Therefore, when g 4 is less 
positive, the g.s. will prefer this phase. Thus the phase 
transitions manifested above reflects the competition of 
the interactions of the three s-channels. The competition 
depends on M because the least number of spins that 
must be aligned is determined by M. 

When M = 2TV is given, the ferro-phase is the only 
choice disregarding how the interaction is as shown by 
the uppermost curve in Fig.l. When TV is sufficiently 
large the terms ~ 1/TV can be neglected, then we have 

E v s ~ (02 + g v + 4 gs)N 2 

= TV 2 [-(cost/? + sinyj) + 1.25lxt(cl>) 1 / 2 . (9) 

6 

From this formula we know that E v g has a minimum 
at p = 7 t /4 and a maximum at p = 57 r /4 where 
E vS = A )N 2 u{Cj) 1 / 2 . Since the total energy is 

proportional to N 3 ^ 5 (E v s) 2 ^ 5 as shown in Eq.( 8 ), with 
TV = 12000 and our parameters, the minimum at p = ir/4 
has E gs = 51612feu, and the maximum at p = 57 r /4 has 
E gs = 60125/ia;. This explains the origin of the dip and 
the highest peak in Fig.l. 

Since the root mean square radius is an observable, 
it is desirable to study R lms so that the theory can be 
checked via experimental measurement. According to 
Eq.(7), .Rrms is proportional to (Egg) 1 / 2 , thus the vari¬ 
ation of i?rms versus p with the minimum and maximum 
is predicted. In particular, when p = 7 r /4 (57 t/4 ) and the 
unit A is replaced by /im, the minimum (maximum) in 
Fig.l is associated with R rms = 3.53 (3.81) pm. In region 
II and III, it is predicted that the increase of M would 
lead to the increase of the size. This is also a point to be 
checked. 


V. STABILITY OF THE GROUND STATE 

It is believed that the stability of the g.s. is assured by 
the gap E gap , namely, the energy difference between the 
first excited state and the g.s.. The former has its (v, S) 
slightly different from that of the g.s., (v g , S g ), and can be 
easily obtained from Eq.(4). An example of A gap is shown 
in Fig.2. Note that when M = 2 TV, (v, S ) = (TV, 2 TV ) is 
the only choice, no choice other than ( TV, 2 TV) is allowed. 
Thus, under M- conservation, there are no excited spin- 
modes. Hence, in order to show the gap for a very large 
M, the curve ”5” in Fig.2 is not given as M — 2 TV but 
M = 2N — 4. Fig.2 has the following features: 

(i) Since Eq.(4) manifests that the energy is nearly 
oc v 2 and S' 2 , the gap would be in general large if both 
v g and S g are large. In this case, a small deviation in 
(v, S ) will cause a great change in energy. Whereas if 
one of them is zero, the gap would be much lower. In 
region I the g.s. is in the ferro-phase with (v g ,S g ) = 
(TV, 2N). Therefore, the gap is very high and accordingly 



FIG. 2: (color on line) The energy gap E gap (in Tko) versus 
i p. The parameters are the same as in Fig.l, except that the 
curve ”5” has M = 2 N — 4 instead of M = 2 N. 


the g.s. is very stable in the ferro-phase. In II and III, 
the upmost curve in Fig.2 has M = 2TV —4. Accordingly, 
(■ v g ,S g ) = (TV — 2, 2TV — 4) and (TV, 2TV — 4), respectively. 
Since they are also large, the gap is also high. When 
M = 2TV — 4, the g.s. is very close to be in ferro-phase 
and therefore is stable. However, when M decreases, the 
stability reduces. In particular, when M = 0, (v g ,S g ) = 
(0, 0) and (TV, 0) in II and III, respectively. In this case 
E ga .p < 0.00012 was found in both II and III as shown 
by the lowest curve in Fig.2. Thus the g.s. in polar and 
cyclic phases with a small M are highly unstable. 

(ii) Egap in every region appears as a peak, namely, in 
the middle part it is higher but very low when p is close to 
the borders. For Fig.2 the g.s. in I has (v, S ) = (TV, 2TV), 
while the first excited state has (v,S) = (TV, 2TV — 2) 
if tp < 0.47227T. Thus the excitation from (TV, 2TV) to 
(TV, 2TV — 2) is caused by a change in S. However, when 
ip —> 0, accordingly gs —> 0. Therefore the gap is zero. 
On the other hand, the first excited state has (v, S) = 
(TV — 2, 2TV — 4) when (0.4722-7T < ip < p a ). From Eq.(4), 
we know the gap 

En,2N — En-2,2N-4 

= ^[7(TV+i) ff0 4 + 10(TV-2) 524 ]. (10) 

One can prove that the gap is zero when p = p a i = 
tan -1 I'u 3 (-jv^- 2 ) • Since p a > is extremely close to p a , this 
explains the decline of the gap when p —> p a . Based 
on Eq.(4), the decline of the peak in II and III can be 
similarly explained. Thus, in the neighborhoods of the 
borders, the g.s. is highly unstable. 

In addition to the gap, another factor that could af¬ 
fect the stability of the g.s. is the level density in the 
neighborhood of the lowest level. This density can be 
calculated based on Eq.( 8 ). As an example, the num¬ 
ber of levels with excitation energy < 0.1 Huj is given in 
Table U in which seven choices of p and four choices of 









M are chosen. The other parameters are the same as in 
Fig.l. If the degeneracy of a level is d , then the number 
contributed by this level is d. 

TABLE I: The number of levels lower than 0.1 hu; versus M 
and p. p a and pb denote the borders, 5 = 0.017T. The first 
column gives the region of p. 


Region 

P 

Number of levels 

M = 0 

5 : 

II 

M = 2N - 4 

I 

7t/4 

1 

1 

1 

I 

Pa ~ S 

4 

4 

2 

II 

Pa + S 

23059 

7 

2 

II 

TV 

20596 

1 

1 

II 

Pb — s 

167171 

19 

2 

III 

Pb + S 

3372 

7999 

2 

III 

3tt/2 

885 

2001 

1 


From the table we know that 

(i) When ip is in region I or M is close to 2 N (i.e., in or 
close to the ferro-phase), the density of low-lying states 
is rather diffuse. Together with the big gap, both factors 
assure the stability of the g.s. (When the number = 1 as 
shown in the second row, there is no excited states lower 
than 0.1 hui). 

(ii) When tp is in region II and M is very small (i.e., 
in or close to the polar phase), the density is extremely 
dense. Together with the very small gap, both factors 
lead to a highly unstable g.s.. This situation can be 
greatly improved when M gets larger. 

(iii) When p is in region III and M is not very large 
(say, M < N), the low-lying spectrum is also dense but 
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not as dense as in the polar phase. Note that, when M is 
not close to 0 or 2TV (say, 0.05 x 2TV < M < 0.95 x 2TV), it 
has been mentioned that the states in cyclic phase could 
be highly degenerate. The great degeneracy contributes 
to the level-density substantially (say, in the last row 
of Tabled the number 2001 arises completely from the 
degeneracy of the g.s.). 


VI. NUMERICAL SOLUTIONS OF THE 
SYMMETRY ADAPTED GP-EQUATION 

By making use of the TFA we have obtained analytical 
solutions, thereby the related physics can be understood 
in an analytical way. In order to evaluate the accuracy of 
the TFA, the equation with the kinetic energy included 
is solved numerically and the results are given below for 
a comparison. Firstly, we found that the curves for E gs 
from the numerical solutions as a whole is an upward 
shift of those plotted in Fig.l. It implies that the amount 
of total kinetic energies contained in various states with 
very different spin-textiles are similar, and the E gs given 
under the TFA are correct when N = 12000 (except the 
omission of the kinetic energy). The shift is about 3.5 x 
I0 3 hoj (thus the kinetic energy is about 6% of the total 
energy in our cases). Secondly, we found that the curves 
for -Egap from the numerical solutions overlap the curves 
from TFA plotted in Fig.2. It implies that the amount 
of kinetic energies contained in the g.s. and in the first 
excited state is similar. Thus the spin-excitation does 
not remarkably affect the spatial motion. 


TABLE II: Root Means Square Radius Arms in pm from nu¬ 
merical calculation versus N, <p, and M with uo = 300s _1 . 
The values from TFA are given inside the parentheses. 


N 

P 

M 

^rms 

N 

P 

M 

Rrms 

N 

P 

M 

•firms 

1200 

7t/4 

0, N, 2 N 

2.48(2.23) 

6000 

7t/4 

0, N, 2 N 

3.19(3.07) 

12000 

7t/4 

0, N, 2 N 

3.62(3.53) 


TV 

0 

2.57(2.35) 


TV 

0 

3.34(3.24) 


TV 

0 

3.80(3.72) 



N 

2.58(2.36) 



N 

3.36(3.25) 



N 

3.81(3.73) 



2 N 

2.60(2.38) 



2 N 

3.39(3.29) 



2 N 

3.85(3.77) 


57t/4 

0 

2.56(2.33) 


57r/4 

0 

3.33(3.22) 


57t/4 

0 

3.78(3.70) 



N 

2.58(2.35) 



N 

3.35(3.25) 



N 

3.80(3.73) 



2 N 

2.62(2.40) 



2 N 

3.42(3.32) 



2 N 

3.88(3.81) 


Furthermore, the size of the system is expected to in¬ 
crease by including the kinetic energy. The results from 


the numerical calculation are shown in Table El where 
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FIG. 3: (color on line) log 10 F? rms versus log 10 N with ui = 
300s _1 where R IU1S is in /.im. M is given at three values, 
i p = 7t/4, it, and 57t/4 for the curves ”1” to ”3”, respectively. 
The lower panels are for the slopes of the curves. 
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FIG. 4: (color on line) log 10 firms versus log 10 Hi with N = 
12000. Refer to Fig.3. Note that the vertical scale in 4d to 4f 
is very small. It implies that the curves in 4a to 4c are very 
close to straight lines. 


the f? rms from TFA are given inside the parentheses and 
are smaller. From the table we know that the TFA leads 
to a ~ 8% reduction of the f? rm s (if AT = 1200), or a 
~ 2% reduction of the fi rms (if N = 12000). Nonethe¬ 
less, the two features found previously under the TFA re¬ 
main unchanged, namely, (i) When the g.s. is not in the 
ferro-phase, a larger M leads to a larger size, (ii) When 
the g.s. is in the ferro-phase (p < p a or M = 2 N), the 


size is minimized when p = it/ 4 and is maximized when 
p = 57t/4. 

It has been predicted based on the TFA that F? rms oc 
N 1 2 3 / 5 . To check this relation log 10 f? rms from the numeri¬ 
cal solutions are plotted in Fig.3 versus log 10 N. In 3a to 
3c the curves are not exactly straight lines. Their slopes 
depend on N and are plotted in 3d to 3f, respectively. 
Disregarding M and p the slopes tend to 1/5 when N 
increases as predicted. It has been predicted based on 
the TFA that F? rms A oc (w) -2 / 5 . To check this relation 
log 10 f? rms from the numerical solutions are plotted in 
Fig.4 versus log 10 w. The slopes of the curves tend to 
—2/5 when Co increases as predicted. 


VII. SUMMARY 


The generalized GP equation adapted to the Cf(5) D 
SO(5) D SO( 3) symmetry has been derived for spin-2 
condensates. This equation has been solved analytically 
under the TFA and by strict numerical calculation. It 
was found that the TFA is applicable if N is large (say, 
N > 10 4 ). The emphasis is placed on the g.s.. Based 
on a rigorous treatment of the spin-degrees of freedom, 
the detailed spin-textiles, i.e., the ferro-, polar, and cyclic 
phases, and their mixing, are explained in a many-body 
way and thus the underlying physics can be understood 
beyond the mean-field-theory. Besides, the variation of 
the spin-textiles versus M in regions II and III is no¬ 
table. Note that the factors affecting the stability of the 
g.s., such as the gap and the neighboring level density, 
together with the degeneracy of the g.s. itself, are less 
touched in existing literatures. These factors are studied 
in detail in this paper. The great difference in the sta¬ 
bility and degeneracy caused by varying p (which marks 
the features of the interaction) and M is notable (this 
is true even when p varies within a region, i.e., the g.s. 
remains in the same phase). We believe that the effect 
of these factors would be serious when the temperature 
is very low. Since f? rms is an observable, efforts have also 
been made to clarify the relation between R rrns and N, 
to, and p. This provides a way for checking the theories 
with experimental data. 
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